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1 INTRODUCTION 



Strong variability of the X-ray flux is a generic property 
of the black hole binaries in the so called hard state (see 
e.g. van der Klis 1994). While it is usually assumed that 
X-ray emission is produced via Comptonization of the low 
frequency photons by an optically thin cloud of hot electrons 
(e.g. Sunyaev and Truemper 1979) the nature of the strong 
variations is less well understood. One of the possible clues 
to the origin of the variability may be connected with the 
presence of the time lags between the light curves in differ- 
ent energy bands. The time lags were first found in the early 
observations of Cygnus X-l (Priedhorsky et al. 1979, Nolan 
et al. 1981) and confirmed in the numerous subsequent ob- 
servations of Cygnus X-l and other black hole candidates. 
In the vast majority of cases hard photons lag with respect 
to the soft ones. From observations it was found that (i) the 
value of the time lag depends on frequency, (ii) time lag ap- 
proximately logarithmically depends on the energy separa- 
tion between two bands (e.g. Miyamoto & Kitamoto 1989), 
(iii) cross-correlation function between two bands peaks at 
zero lags within ~1 ms (Maccarone et al. 2000 and refer- 
ences therein). Comprehensive description of the observed 
Cygnus X-l variability and time lags are given by Nowak et 
al., 1999a,b and Pottschmidt et al. 2000. A recent reviews 



ABSTRACT 

It is shown that the energy dependence of the time lags in Cygnus X-l excludes 
any significant contribution of the standard reflected component to the observed lags. 
The conclusion is valid in the 0.1-10 Hz frequency range where time lags have been 
detected with high enough significance. In fact the data hint that reflected component 
is working in opposite direction, reducing the lags at energies where contribution of 
the reflected component is significant. 

We argue that the observed logarithmic dependence of time lags on energy can 
be understood as due to the small variations of the spectrum power law index in a 
very simple phcnomcnological model. We assume that an optically thin flow/corona, 
emitting a power law like spectrum, is present at a range of distances from the compact 
object. The slope of the locally emitted spectrum is a function of distance with the 
hardest spectrum emitted in the innermost region. If perturbations with different time 
scales are introduced to the accretion flow at different radii the observed X-ray lags 
naturally appear due to the inward propagation of perturbations on the diffusion time 
scales. 
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of the observed lags and theoretical models are given by 
Poutanen 2001a,b. 

We discuss below the time lags using RXTE data on 
Cygnus X-l as the example. In Section 2 we show that ob- 
served energy dependence of the time lags excludes signifi- 
cant contribution of the standard reflected component to the 
time lags in Cygnus X-l. In Section 3 we introduce simple 
phenomenological model of the time lags, which naturally 
explains the logarithmic energy dependence of the time lags. 
In Section 4 we discuss the possibility of explaining the lags 
as due to diffusive propagation of the perturbations in the 
accretion flow. In Section 5 we discuss implications of the 
results. Section 6 summarizes our finding. 



2 TIME LAGS DUE TO LIGHT CROSSING 
TIME OF THE REFLECTOR 

For hard X-ray state of Cygnus X-l one of the popular 
models assumes that an optically thick and geometrically 
thin disk is truncated at some radius from the compact 
object and most of the observed X-ray emission is com- 
ing from the more compact optically thin region located 
closer to the black hole (e.g. Esin, McClintock and Narayan 
1997, Zdziarski et al. 1999). A fraction of the hard X-ray 
photons emitted by the central region can be intercepted 
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by the optically thick disk and reemitted again. This re- 
emitted/reprocessed component, also called "reflected" com- 
ponent, has a very distinct spectral shape (e.g. Basko, Sun- 
yaev & Titarchuk 1974, George & Fabian 1991). If matter 
in the optically thick disk is neutral or weakly ionized than 
the spectrum of the reflected component is harder (in the 
standard X-ray band) than the illuminating spectrum and 
it also contains prominent iron fluorescent K a line. Because 
of the delay of photons which traveled from the emission re- 
gion to the reflector and then to the observer hard lags may 
naturally appear (e.g. Poutanen 2001b). 

In the presence of the direct and reflected components 
the observed spectrum S(E, t) at time t can be written as: 



S(E,t) = D(E,t) + R(E,t) 



(1) 



where D(E, t) and R(E, t) are the time dependent spectra of 
the direct and reflected components respectively. Assuming 
that the direct component vary in intensity only and not in 
shape: 



D(E,t) =A(t)D(E) 



(2) 



and neglecting the dependence of the reflected component 
spectral shape on the incident and emission angles the re- 
flected component can be written as: 



R(E,t) = ^R(E) 



T(r)A(t - r)dT 



(3) 



where R(E) are the energy spectra of the reflected com- 
ponent calculated for the reflector subtending solid angle 
of 2-7T, A(t) is the normalization of the direct component 
and T(t) is the reflected component transfer function, which 
characterizes the time dependence of the intensity of the re- 
flected component in response for a infinitely short flare of 
the direct component. In the above expression T(t) is nor- 
malized such that J T(t)dt = 1 and the contribution of 
the reflected component to the time averaged spectrum is 
expressed through the standard factor where f2 is the 
total solid angle subtended by the reflector. 
The expression (1) can be rewritten as: 



S(E,t) = A(t)D(E) + £-R(E) 

Z7T 



T(j)A(t - r)dr 



(4) 



Further assuming that contribution of the reflected 
component to the total spectrum is small one can write the 
expression for a frequency resolved spectrum, introduced in 
Revnivtsev et al. 1999: 



\A(f)\D(E) x 



S(E,f) = \S(E,f)\ 
2tt D(E) 1 yJn 



\A(f)\ x (D(E) + ^-R(E)Re[T(f)]) 



(5) 



The equivalent width of the iron line in the frequency re- 
solved spectrum is then obviously: 



EW(f)nEWtxRe[T(f)], 



(6) 



where EWt is the equivalent width in the total source spec- 
trum. From the above expression it is clear that the equiv- 
alent width in the frequency resolved spectrum should de- 
crease if the transfer function is smooth at a given frequency. 
E.g. when the light crossing time of the reflector is large the 
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Figure 1. The energy dependence of the time lag in Cygnus X-l 
in the hard state at the frequency of 2.5 Hz. For comparison the 
energy dependence of the time lags expected in the model with 
extended reflector is shown. For upper curve we used relative 
normalization of the Gaussian and reflected continuum obtained 
from the fit to the total Cygnus X-l spectrum during observation 
10238-01-08-00 (Gilfanov et al., 1999). For the lower curve the 
normalization of the Gaussian was chosen such that equivalent 
width of the line with respect to the reflected continuum is ~1 
keV. The latter value is characteristic for an angle averaged re- 
flection from a neutral matter with a normal abundance of heavy 
elements. 



equivalent width should be low at the corresponding time 
scales. Using these arguments Revnivtsev et al. 1999, Gil- 
fanov et al. 2000 set an upper limit on the effective size of 
the reflector in Cygnus X-l of ~50 R g . We note here that 
the above expression assumes linear relation between vari- 
abilities in the reflected and direct components in the form 
of the convolution (^). If on the contrary the variability of 
the direct and reflected components are completely incoher- 
ent/uncorrelated then the contribution of the reflected com- 
ponent to the frequency resolved spectrum is proportional 
to square of the solid angle subtended by the reflector (J^r) 2 . 

One can further test the hypothesis of the large effective 
size of the reflector using the time lags. Assuming again that 
the light curve in any energy band can be expressed through 
(^) the expected phase lag (j>(E, f) between given energy E 
and much lower energy (where contribution of the reflected 
component is negligible and the spectrum is solely due to 
direct component) can be written as: 



tg(4>(E,f)) = 
&R(E)Im[f\ 



Im[S*(E, f)D(EJ)} 
Re[S*(EJ)D(E,f)] 



n r(e) 



Im[f(f)] 



(7) 



D(E) + £R(E)Re[f] 2tt D(E) 

Therefore in this approximation the frequency depen- 
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Figure 2. Transfer functions ^-T(t) for an isotropic point source 
of primary radiation located at h = lOR g above a flat or concave 
disks with inner an radius of = 10 and 50-f? 9 and inclination 
of 0°,50°. 

dence of the phase/time lags is due to imaginary part of 
the Fourier transform of the transfer function. Energy de- 
pendence of the lags is simply the ratio of the reflected and 
direct spectral components. Since D(E) is assumed to be 
smooth function of energy the time/phase lag as the func- 
tion of energy should contain prominent iron line with the 
equivalent width comparable with the equivalent width of 
the line in the reflected component, i.e. much larger equiva- 
lent width than in the total spectrum. 

In order to assess possible contribution of the reflected 
component in the observed time lags in Cygnus X-l we cal- 
culated lags in the narrow energy channels. For our anal- 
ysis we use publicly available data of Rossi X-ray Tim- 
ing Explorer observations P10238 performed between Mar. 
26, 1996 and Mar. 31, 1996 with total exposure time ~70 
ksec (we used only observations with all 5 PCU turned on). 
We used PCA data in the "Generic Binned" mode, having 
jj^sec~16 msec time resolution in 64 energy channels cover- 
ing the whole PCA energy band (B_16ms_64M_0_249). The 
lags were calculated for every channel above 3.7 keV, with 
respect to the count rate in the 2.8-3.7 keV energy band. The 
resulting energy dependence of time lag at the frequencies 
of the order of 2.5 Hz is shown in FigJin. For comparison we 
plot in the same figure an expected energy dependence of 
the lags due to the reflected component, according to equa- 
tion (T7j). Normalization of the expected curve is arbitrary. 
Here we used as D(E) the averaged spectrum of the source 
during this observation. Reflected component R(E) was cal- 

* Time lags at other frequencies in the range from ~0.6 to ~10 
Hz have approximately similar energy dependence (see Fig.H) 



culated using XSPEC Vll (Arnaud 1996) model pexrav 
of Magdziarz & Zdziarski 1995 plus a Gaussian at the en- 
ergy of 6.4 keV. The resulting spectrum has been convolved 
with the PCA energy response matrix. Note that expected 
lags were calculated using equation (Q), i.e. with respect to 
the energy range where contribution of the reflected com- 
ponent was assumed to be zero. Therefore expected lags do 
not go to zero at the energy <4 keV. Two curves shown in 
Fig.nl differ in the normalization of the Gaussian added to 
the reflected continuum. For upper curve we used relative 
normalization of the Gaussian and reflected continuum ob- 
tained from the fit to the total Cygnus X-l spectrum during 
observation 10238-01-08-00 (Gilfanov et al., 1999). For the 
lower curve the normalization of the Gaussian was chosen 
such that the equivalent width of the line with respect to the 
reflected continuum is ~1 keV. The latter value is character- 
istic for an angle averaged reflection from a neutral matter 
with a normal abundance of heavy elements. 

One can see from Fig.|l] that the lags due to the reflected 
component should contain a prominent feature at the energy 
~6.4 keV. We stress again that this feature should be much 
more prominent in the lags than in the averaged spectrum, 
where iron line provides at most 10% relative deviation of 
the observed spectrum from a power law fit. Clearly no such 
feature is present in the observed lags. We therefore conclude 
that the observed lags are predominantly not due to the 
reflected component. 

One can use the same argument in order to constrain 
geometry of the reflecting region. Indeed lack of the iron line 
feature in the observed energy dependence of the time lags 
imply certain constrains on the combination of the reflected 
component strength and the duration of the time delay due 
to the finite size of the reflector. In the Fig.^ we show the 
observed energy dependence of the time lags at three repre- 
senting frequencies: 0.6, 2.5 and 10 Hz. For comparison we 
show in the same plot expected time lags for different geome- 
tries of the reflector and different strength of the reflected 
component. Four basic models have been considered: 

• A: A point isotropic source above the flat disk with the 
hole in the middle. The height h of the source above the disk 
plane is set to 10R g , the size of the hole Ri n in the disk is 
set to 50-R 9 . Here R g is the gravitational radius for a WMq 
black hole. The factor Q./2-k, calculated for this geometry, is 
~ 0.2. 

• B: The same geometry as for A, but the strength of the 
reflected component was enhanced, so that effective param- 
eter Q./2-k is 0.54. This is possible if e.g. the source of the 
primary emission is anisotropic and more emission is coming 
towards the disk than to the observer. Particular value of the 
normalization of the reflected component in this model was 
taken from the fit to the total Cygnus X-l spectrum during 
observation 10238-01-08-00 (Gilfanov et al., 1999). 

• C: The same geometry as for A, but for Ri n = 
10i? 9 .The factor Q/2tt, calculated for this geometry, is ~ 
0.7 

The transfer functions in the form ^-T(t) for models A 
r-i 2lr v ' 

and C are shown in FigI2J for two values of the disk inclina- 
tion 9 (0° and 50°). Model B obviously has the same shape of 
the transfer function as model A and is different only in nor- 
malization. For comparison we also show the transfer func- 
tion for a concave disk with the dependence of disk height 
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Figure 3. The energy dependence of the time lag in Cygnus X- 
1 in the hard state at different frequencies. For comparison we 
show the expected curves for the models where reflected compo- 
nent originates from a flat disk with a central hole illuminated 
by the compact source located at the height h = !0R g above the 
disk surface. The model curves were calculated for the equivalent 
width of the Fe K a line obtained from the fit of the time averaged 
source spectrum. 



over radius adopted from Shakura and Sunyaev (1973). For 
the frequency range considered here such concave disks do 
not differ from the fiat disks. 

The lags expected in these models (for inclination 9 = 
50°) are shown in Fig.|§| The model curves were calculated 
for the equivalent width of the Fe K a line with respect to 
the reflected continuum obtained from the fit of the time 
averaged source spectrum. As was mentioned above lack 
of prominent iron feature in the energy dependence of the 
time lags suggests that lags are largely not due to the re- 
flection. Moreover lags expected in models B and C at 10 
Hz must produce clear feature at 6.4 keV in contradiction 



to the data. Even for model A expected lags can be only 
marginally "hidden" in the observed lags. There are how- 
ever many uncertainties and simplifications in the models 
discussed above. In particular for simplicity we neglected 
the dependence of the reflected component spectrum on the 
incident angle. Also the surface of the reflector (upper few 
Thomson optical depths) may not be smooth and flat, but 
wrinkled. Furthermore ionization of the reflector surface will 
modify the spectrum of the reflected component. Finally all 
relativistic effects associated with reflection from the inner 
regions of the accretion disk (see e.g. Campana and Stella 
1995) have been ignored. Account for all these effects may 
modify the normalization and shape of the energy depen- 
dence of the time lags by a factor of few. 

Perhaps the most dubious assumption is the "linear" 
response of the reflected component to the variations of 
the illuminating flux. Recent analysis (e.g. Nayakshin and 
Kallman 2000) shows that shape of the reflected component 
can vary in a complicated fashion in response to increase of 
the primary flux. Therefore conservative conclusions of this 
section are that i) "linear" reflected component, containing 
prominent fluorescent line at 6.4 keV, is not responsible for 
the observed lags and ii) some of the source/reflector geo- 
metrical models can be definitely excluded. 



3 TIME LAGS DUE TO VARIATIONS OF THE 
POWER LAW INDEX 

The spectrum of Cygnus X-l in the 2-20 keV energy band 
can be reasonably well approximated by a power law. De- 
viations, in particular due to the reflected component, are 
present, but their amplitude relative to the power law com- 
ponent is usually less than 10% in this energy range. We 
assume below that at any moment of time the spectrum can 
be represented as a power law: 



S(E,t) = A(t)E 



-«(*) 



(8) 



where A(t) is the normalization of the power law as a func- 
tion of time and a(t) is the photon index, which is also a 
function time. For simplicity we set a(t) = ao — 0(i), where 
ao is a photon index of the time averaged spectrum and j3(t) 
is the time variable part of the photon index. Assuming that 
/3(t) <S ao the spectrum of the source can be rewritten as: 



S(E, t) = A(t)E~ a ° (1 + 0(t)lnE) 



(9) 



The phase lag cf> at a given frequency / can then be writ- 
ten through the Fourier transforms S* (E, f) of the source 
light curves at two different energies E\ and E2 as: 

. ^ Jm[S*(gi,/)g(g a ,/)] 
Im \{A* + (PA)* lnE 1 )(A+ {p'A)lnE 2 ) 



A* A 



Im 



A' [13 A) 



A* A 



, E 2 
Xln El 



(10) 



In the above expression we neglected terms of the order of 
0(/3 2 ). Thus the assumption that the spectrum has a power 
law shape at any moment of time and the variations of slope 
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Figure 4. Time lag at the frequency of 2.5Hz as a function of 
Log(E). See text for the description of the solid line. 



are small automatically implies that phase lags (and conse- 
quently time lags) scales as In^-. 

3.1 Universal slope flux correlation? 

It is clear that if variations of slope j3(t) and normalization 
A(t) are random and incoherent then the time lag is zero. 
If these two quantities are correlated then nonzero lag may 
appear. The simplest assumption would be a linear corre- 
lation between slope and normalization (i.e. between slope 
and flux): 



P{t) = 7 x {A{t) - A ) 



(11) 



where Ao is the normalization of the time averaged spectrum 
and 7 is the coefficient of linear proportionality between 
variations of normalization and slope. 

The expression for phase lags is then further reduces to: 



Im A* (A 2 ) 



A* A 



, E 2 
x "/In — 
E\ 



(12) 



In this approximation the sign and the amplitude of the 
lags would simply reflect properties of the flux variations 
on different time scales. Unfortunately the function (|l2|), 
calculated for the observed light curves, changes the sign 
several times within the frequency range where hard lags are 
observed. This means that simple linear correlation between 
the flux and the slope of the spectrum can be excluded or 
that sign of the correlations between flux and slope also has 
to change sign in order to keep the sign of ( |l2"| ) the same at 
all frequencies. We note here that change of the slope/flux 
correlations sign at different time scales was indeed reported 
for Cygnus X-l (Li, Feng & Chen 1999). Anyway it is clear 
that simple linear correlation between the slope and the flux 



fails to produce observed hard lags and more complicated 
models are to be invoked. 



4 TIME LAGS DUE TO ACCRETION 

In this section we construct simple model which attributes 
the lags to the propagation of perturbations in the accre- 
tion flow from the outer to inner region. As was discussed 
by Lyubarskii (1997) and Churazov, Gilfanov and Revnivt- 
sev (2001) strong variability of the X-ray flux over very 
broad dynamic range of time scales suggests that different 
time scales are introduced to the accretion flow at differ- 
ent distances from the compact source. One might think 
e.g. in terms of the effect on the mass accretion rate of the 
magneto- hydrodynamic turbulence, which serves as a source 
of the viscosity in the accretion flow through the fluctuating 
magnetic stresses (e.g. Balbus and Hawley 1991, Hawley, 
Gammie and Balbus 1995, Brandenburg, Nordlund, Stein 
and Torkelsson 1995). The perturbations introduced at large 
distance from the compact object are then propagated down 
to the region of main energy release and cause observed vari- 
ations of the X-ray flux. If the propagation of the perturba- 
tions in radial direction is due to viscous diffusion then the 
very fact that we observe strong variations of the X-ray flux 
means that propagation/diffusion time from a given radius 
ro is comparable to or shorter than the time scale of pertur- 
bations, introduced to the flow at this radius ro. Otherwise 
these perturbations would be completely washed out before 
reaching the innermost region. This in turn means that lags 
could naturally appear in such situation if locally emitted 
spectrum is a function of radius r. Of course energy budget 
of the accretion flow far from the innermost region is small, 
but the observed lags are also small, of the order of 0.1 ra- 
dian. We then tested feasibility of this model using following 
assumptions: 

• Statistically independent (incoherent) perturbations 
are introduced to the accretion flow at different distances 
from the compact source. The characteristic time scale of 
these perturbations is a function of radius. 

• We further assume that the shape of the locally emitted 
spectra depends on radius and it is getting progressively 
softer with increase of the distance from the compact object. 



The light curve at a given energy E is: 



L(E,t) 



D(r ,t )G(r, r ,t- t )e(r)S(E, r)drdr dt (13) 



where D(ro,to) is an initial perturbation introduced to the 
flow at radius ro and time to, G(r,ro,t — to) is the Green 
function, which describes propagation of the perturbations 
from radius ro to r, S(E, r) is the shape of the locally emit- 
ted spectrum at the radius r and e(r) is the total luminosity 
of the accretion flow emitted at a given radius. Assumption 
of statistical independence of perturbations introduced at 
different radii implies that quantities like power density or 
cross spectra can be calculated independently for each ini- 
tial radius and then averaged over range of initial radii. The 
Fourier transform L(E, /) of the X-ray light curve is then: 



L(E,f) 



D(r ,f) / G(r,r ,f)e(r)S(E,r)drdr (14) 
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where D(ro, /) is the Fourier transform of the perturbations 
introduced to the flow at the radius ro, G(r,ro,f) is the 
Fourier transform of the Green function, / is the frequency. 
From the above equation it is clear that in this approxima- 
tion time lags are solely due to the properties of the Green 
function and the dependence of the spectral shape on radius. 
For our toy model we choose the simplest possible version 
with e(r) oc 1/r 2 per unit dr. In order to calculate phase 
and time lags we created two artificial light curves in the 
"soft" and "hard" bands assuming simple laws of emissivity 
in these bands: 



S(r)=e(r)S(E soft ,r) = - 

H(r) = e(r)S(E hard ,r) = \x h(r) 



(15) 



where h(r) is a decreasing function of the radius. For sim- 
plicity we assume that propagation of the perturbations can 
be characterized by the Green function of Lynden-Bell and 
Pringle (1974) describing diffusion in the geometrically thin 
disk approximation (see also Lyubarskii 1997). I.e. 



G(r,r ,t) oc 



-1/4 



-exp 



X/IX 



+ r 



1/21 



it 



-l/4i 1/4! 



■It 



(16) 



where I v (x) is the Bessel function of imaginary argument. 
For the parameter I we use the value of 1/3 appropriate 
for a standard a disk with constant ratio — , where H is 
a half thickness of the disk. In this equation and below we 
measure radius r in units of r m in, which is supposed to be 
the characteristic radius of the inner region of the accretion 
flow where most of the energy is released. Time is measured 
in units of the characteristic diffusion time k 2 x r^J 2 ^ , where 
k is the effective diffusion coefficient. 

Two simple versions of the function h(r) (hereafter 
"hardness" function) are shown in Fig.^. The dotted line 
shows the case when slope of the locally emitted spectrum 
is the same at all distances from the compact object, but it 
abruptly hardens within innermost region. The solid line 
shows the case when locally emitted spectrum gradually 
hardens as the distance from the compact object decreases. 
In Figj^ light curves in the soft and hard energy bands are 
shown, assuming h(r) function with the abrupt jump. Here 
initial radius ro was set to 30 and initial perturbation was 
assumed to be S function in time and radius. Initially lo- 
cally emitted spectrum is soft and it hardens when pertur- 
bation spreads all the way down to the minimal radius. After 
that moment hardness of spectrum remains practically un- 
changed. 

We then calculated phase lags fixing the initial radius 
ro = 30 and varying the characteristic frequency of the ini- 
tial perturbation. The dependence of the phase lag on fre- 
quency is shown in the Fig.j?]. Solid line in this figure shows 
the phase lag (in radians) as a function of frequency ex- 
pressed in units of inverse characteristic diffusion time for a 
given initial radius r^ 2 . These lags were calculated for h(r) 
function with a jump (see Figjjs]), but for the smoother h(r) 
the resulting lags are practically the same. Dotted line in the 
same figure shows the value of the transfer function power 
| J G(r, ro, f)e(r)dr\ 2 . It characterizes to what extent an am- 
plitude of a perturbation at a given frequency / is suppressed 
in the observed light curve due to diffusive spreading. Val- 
ues of the order of 10 -2 imply that "local" perturbations 




Figure 5. Two versions of h(r) function, used in the simulations. 
Dotted line shows the h(r) function when the locally emitted 
spectrum has the same slope at all distances from the compact 
object except for the innermost region, where the spectrum is 
harder. Solid line shows the case when spectrum gradually steep- 
ens with distance from the compact object. The amplitude of the 
h(r) function change affects the absolute value of the time/phase 




Figure 6. Light curves in two energy bands (soft — solid line, 
hard — dashed line) for the <5(t) perturbation at the initial radius 
r = 30. 
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Figure 7. Phase lag (solid line) as a function frequency for 
the perturbations introduced at fixed initial radius ro = 30. 
Dotted line show the value of the transfer function power 
| J G(r, ro, f)e(r)dr\ 2 , which characterizes the suppression of the 
power of the initial perturbation at a given frequency, due to dif- 
fusive spreading. 



Figure 8. Phase lags for two versions of h(r). Frequency, intro- 
duced at radius ro is assumed to scale as K X r. 



-3/2 



with K = 1. 



at ro with an amplitude of the order of unity would result 
in ~10% variations in the innermost region where bulk of 
energy is released. 

In order to calculate phase lags associated with the per- 
turbations generated at different radii we made further sim- 
plifying assumption that at every initial radius ro instabili- 
ties generate single frequency which scales as a characteristic 
diffusion time at this radiusQ I.e. f(ro) = K * r g 3 ^ 2 , where 
if is a parameter. From Fig.|^ one would expect the lags 
to increase with the increase of K unless K is too large. 
The resulting phase lags for two versions of the "hardness" 
function are shown in the Fig.^. The resulting phase lags 
weakly grows with frequency and reach ~ 0.1 radian before 
falling sharply. The amplitude of the phase lags and the de- 
pendence of the lags on frequency are thus regulated by two 
major parameters: i) the strength of the spectral changes as 
a function of radius and ii) the relation between the time 
scale of perturbations, locally generated at radius ro, and 
diffusion time scales from ro down to much smaller radii. 
We stress that under assumption that at each radius ro per- 
turbations with single characteristic frequency are generated 
the phase lags do not depend on the shape of the power den- 
sity spectra and are solely determined by the Green function 
and the "hardness" function. 




0.001 



0.01 0.1 
Frequency, f=K r^ 3/z 



le do 
3© 



soft light curves (see eqJla) for two versions of h(r). 



t Note that for the a disk with constant ratio of the disk height 
to the radius the characteristic diffusion time and Keplerian time 
scales similarly 
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Figure 10. Lags for the free fall model for K = 1 and K = 0.5. 

In the model, considered above, one would expect that 
power density spectra in different energy bands should have 
somewhat different shape due to the influence of the "hard- 
ness" function. In the Fig|^ we show the ratio of the power 
density spectra in the "hard" and "soft" bands as a func- 
tion of frequency. This ratio is again independent on the 
shape of the power density spectrum and is governed by the 
the Green function and the "hardness" function. One can 
see that, as expected, the ratio of the power density spectra 
increases with frequency. 

For completeness we also calculated time lags using the 
same choice of the "hardness" function, but instead of dif- 
fusion Green function we substitute Green function appro- 
priate for a free fall from given radius ro. This is equiva- 
lent to the assumption that distinct blobs of emitting mat- 
ter are falling towards the compact source with the velocity 
v oc r- _1//2 . The resulting lags are shown in Fig.|lo|. Since in 
this case no viscous spreading is present, the absolute value 
of the lags is higher for the same choice of the "hardness" 
function and parameter K. 

Finally we stress that if the locally emitted spectrum 
at any radii is a power law then observed lags will show 
logarithmic dependence on energy even if the total observed 
spectrum (integrated over all radii) is not well represented 
by a power law as assumed in the previous section. This 
follows from our assumption that perturbations generated 
at different radii are statistically independent. 

4.1 "Antilags" due to reflected component 

As we see in Section 2 observed energy dependence of the 
time lags contradicts the assumption that hard lags are pre- 
dominantly produced by the time delay of the reflected com- 
ponent with respect to the primary continuum. In partic- 



ular (see Figj4|) instead of a prominent hump at the en- 
ergy of the iron fluorescent line the observed lags seem to 
be slightly suppressed here. This behavior might hint that 
reflected component somehow works in opposite direction 
and reduces the lags. We now discuss how such suppression 
could be introduced within the frame of the model consid- 
ered above. 

An assumption that locally emitted spectrum has a 
power law shape with an index slowly varying with the dis- 
tance (so that hardest spectra are near the center) from the 
compact object will naturally lead to a hard lags with a log- 
arithmic dependence of lags on energy. Let us now assume 
that locally emitted spectra also contain a reflected compo- 
nent. Assuming that overall geometry of the accretion flow 
resembles a truncated optically thick disk followed by an 
optically thin flow one would expect that the reflection frac- 
tion is higher for spectra emitted at larger distances from 
the compact object. I.e. softer spectra (emitted at larger dis- 
tance from the compact object) have larger reflection frac- 
tion than harder spectra. In our model this means that re- 
flected component is leading the harder spectrum. In this 
situation one would expect slight modification of the loga- 
rithmic dependence of lags on energy due to reflected compo- 
nent. The expected shape of the lags dependence on energy 
in this case is shown in Fig.ji] with the solid line. For simplic- 
ity we assumed that lags are due to two power law spectra 
with the power law indices of 2.4 and 1.8 and the reflection 
fraction of ~1 and ~0.5 respectively. I.e. the light curve is: 

L(E,t) = Si(E) *S(t) +S 2 (E) *S(t- At) (17) 

Reflected components in both spectra were smeared with a 
~0.8 keV Gaussian. In the limit of small frequencies one can 
easily write an expression for the time lag as a function of 
energy for the light curve described by eq. (|l^). Detailed 
shape of the energy dependence is affected by strength of 
the reflection features in both spectra, amount of smearing 
applied, relative normalization of both spectra. One illustra- 
tive example, calculated for one particular choice of these 
parameters, is shown in Fig.^J with a solid curve. This curve 
is not supposed to closely reproduce the data, but just to 
demonstrate that the trend is in right direction. 



5 DISCUSSION 

The properties of the hard lags, observed in Cyg X-l can be 
summaries as following: 

• In the hard state of Cygnus X-l lags are observed very 
significantly in the frequency range of few 10" 2 -few 10 Hz. 
Phase lag is a slightly growing function of frequency with a 
maximum phase lag of 0.1 rad between 3 and 10 keV (see 
recent data of Nowak et al. 1999a, Pottschmidt et al. 2000, 
also Fig|ll]). 

• An energy dependence of the lags is approximately log- 
arithmic (Miyamoto et al. 1988, Nowak et al. 1999a, see also 
Fig.^) with the slight, but statistically significant, deviations 
from this dependence. 

• The cross correlation function peaks at zero lag (Mac- 
carone et al. 2000). 

• The coherence function is close to unity (Nowak et al. 
1999a), at least in the middle of the frequency range, where 
lags are observed with high significance 
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Figure 11. Cyg X-l: observed power density spectrum (lower 
panel), the phase lags (middle panel) and ratio of the power den- 
sity spectra in the two energy bands (upper panel). 



Several different branches of models have been suggested 
as an explanation for the observed hard lags (see Poutanen 
2001a,b for recent reviews). Our model belongs to the class 
of "propagation" models as defined by Nowak et al. 1999b. 
Other models involving "propagation" are those of Kato 
1989, Nowak et al. 1999b, Bottcher and Liang 1999, Misra 
2000, (see also Manmoto et al., 1996). In fact our version 
of the model is close to the original suggestion of Miyamoto 
et al. 1988, Miyamoto and Kitamoto 1989 that a clump of 
matter ... drifts from the outer (soft-X-ray-emitting) region 
to the inner (hard-X-ray emitting) region. 

The model is based on two major assumptions: a) dif- 
ferent time scales are introduced to the accretion flow at dif- 
ferent distances from the compact object and propagation 
time scales are comparable to the time scales of perturba- 
tions and b) locally emitted spectra are power laws with the 
photon indices smaller in the innermost region. 

The logarithmic energy dependence of the lags follows 
directly from the assumption (b). Zero shift of the cross 
correlation function and coherence function close to unity 
naturally appear because most of the emission is coming 



from the innermost region and therefore most of the flux 
at all energies is released at the same time (see Fig^). The 
dependence of the lags on frequency and maximal values 
of the lags are tunable parameters of the model (through 
the dependence of the photon index on the radius - the 
"hardness" function h(r) and through the relation of the 
perturbation and propagation time scales) . Examples shown 
above (see e.g. Figj^) show that with reasonable assumptions 
one can approximately reproduce both the dependence of 
lags on frequency and the maximal value of the phase lag. 

In addition this model provides natural explanation of 
the ratio of the power density spectra in different energy 
bands as a function of frequency (see Fig.^ and [ll], upper 
panel) and offers a qualitative explanation for the deviations 
of the energy dependence of the lags from a simple logarith- 
mic law (Fig.^J) . Broadly the model is consistent with the 
geometrical model in which the inside truncated optically 
thick disk is followed by the optically thin flow. At the dis- 
tances larger than the disk truncation radius an optically 
thin flow seems to be also present and perhaps has a form 
of a corona sandwiching the disk (see arguments in Chura- 
zov et al. 2001). The latter assumption allows one to extend 
this model to the soft state of Cyg X-l without any modifi- 
cations. 

Our assumption that each frequency in the observed 
flux variability is associated with some particular radius is 
of course a gross oversimplification. However due to sup- 
pression of the high frequency variations in the course of 
diffusion such situation may in practice materialize in the 
accretion flow if the times scales of locally induced pertur- 
bations at given radius are high, e.g. higher than the dif- 
fusion frequency (Churazov et al. 2001). From Fig.jn] it is 
clear that there are several "humps" in the power density 
spectrum of Cyg X-l. Assuming that this is an indication 
that perturbations arising at some particular radius of the 
accretion flow are dominating at a range of frequencies one 
would then expect the phase lags to rise faster with fre- 
quency (compare Fig.[] and Q) . It seems that such behavior 
is indeed observed as a distinct regions of fast phase lag rise 
(or time lag shoulders - see Nowak et al. 1999a,b) in the 
Cyg X-l data (Fig©. 

Particular emission mechanism is not very important as 
long as only the origin of lags is considered. E.g. emission 
could be due to magnetic flares as in the picture of Galeev, 
Rosner, Viana 1979 or Poutanen and Fabian 1999. In our 
model flares spectra need not to evolve (from soft to hard) 
themselves. They could be short and could occur at different 
distances from the compact object. The observed variability 
could then be not due to individual flares but rather due 
to change of the number of active flares at a given moment 
of time. The model requires however that the spectra of a 
more distant flares to be softer than those of the flares in 
the innermost region of the accretion flow. 



6 CONCLUSIONS 

We show that some properties of the black hole candidates 
spectral variability in the X-ray band (in particular time 
lags) can be explained by a simple phenomenological model. 
This model assumes that a) different time scales are intro- 
duced to the accretion flow at different distances from the 
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compact object and propagation time scales are comparable 
to the time scales of perturbations and b) locally emitted 
spectra are power laws with the photon indices smaller in 
the innermost region. 

We also show that energy dependence of the time lags 
excludes the possibility that lags (in the 0.1-10 Hz range) 
are predominantly due to the delay of the reflected compo- 
nent with respect to the illuminating continuum. This con- 
clusion is derived under assumption of a "standard" neutral 
reflected component which linearly responds to the varia- 
tions of the illuminating flux. 
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